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Abstract 

Extremely anisotropic metal-dielectric multilayer metamaterials are designed to have the effective 
permittivity tensor of a transverse component (parallel to the interfaces of the multilayer) with 
zero real part and a longitudinal component (normal to the interfaces of the multilayer) with ultra- 
large imaginary part at the same wavelength, including the optical nonlocality analysis based on the 
transfer-matrix method. The diffraction-free deep-subwavelength optical beam propagation with 
near-zero phase variation in the designed multilayer stack due to the near-flat iso-frequency contour 
is demonstrated and analyzed, including the effects of the multilayer period and the material loss. 
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I. INTRODUCTION 


Light beam propagating beyond diffraction limit with deep-snbwavelength conhnement 
and near-zero phase variation is highly desirable in many optical integration applications. 
In order to realize snch intriguing optical beam propagation phenomena, optical materials 
are required to possess an extremely anisotropic effective refractive index tensor with an 
ultra-high component perpendicular to the propagation direction to support the light beam 
of large wave vector and an ultra-low component in the propagation direction to reduce 
the phase variation of the light beam simultaneously. Optical metamaterials with artihcial 
subwavelength meta-atoms and tunable effective permittivity and effective permeability have 
been designed to exhibit ultra-high effective refractive index so as to achieve large wave 
vectors [1, 2], as well as to have zero effective refractive index for obtaining the optical 
beam propagation with zero phase variation [3, 4]. Among all these designs, the metal- 
dielectric multilayer metamaterials are of great interest in demonstrating many exciting 
applications such as negative refraction [5], epsilon-near-zero (ENZ) materials [6-9], and 
epsilon-near-pole (ENP) materials [10]. Particularly, possessing an indehnite permittivity 
tensor, the multilayer metamaterials are of unique applications in subwavelength imaging 
[11, 12], enhanced photonic density of states [13], broadband light absorbers [14-16], ultra- 
high refractive indices for subwavelength optical waveguides [17], indehnite cavities [18, 19], 
broadband spontaneous emission engineering [20-22], super-Plank thermal emission [23-25], 
ultra-large-wavevector plasmon polaritons [26-28], and other exceptional properties [29-31]. 

In our previous work [32] , the extremely loss-anisotropic metal-dielectric multilayer meta¬ 
materials are designed in order to realize the diffraction-free optical beam propagation mainly 
based on the effective medium theory (EMT) and the operation wavelength is adjusted from 
the dispersion relation of the multilayer stack. However, the effective permittivity ten¬ 
sor including the optical nonlocality is not well considered in the analysis. Therefore, in 
the present work, a new kind of metal-dielectric multilayer metamaterials are designed to 
have the extremely anisotropic effective permittivity tensor for the demonstration of both 
the diffraction-free optical beam propagation and the near-zero phase variation along the 
propagation direction, based on a theoretical analysis of the effective permittivity tensor 
including the optical nonlocality. The designed multilayer metamaterials will have the ef¬ 
fective permittivity tensor with a near-zero transverse component parallel to the multilayer 
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and an nltra-large longitndinal component in the imaginary part normal to the mnltilayer 
at the same wavelength, which represents the ENZ-ENP wavelength. The near-zero trans¬ 
verse component of the effective permittivity tensor leads to the near-zero phase variation 
dnring the propagation of the electromagnetic wave. On the other hand, the longitndinal 
component with an nltra-large imaginary part of the effective permittivity tensor will re- 
dnce the diffraction of the electromagnetic wave [33, 34]. The hlling ratio of metal in the 
metal-dielectric mnltilayer stack with respect to the ENZ-ENP wavelength is fnlly analyzed 
according to the transfer-matrix method inclnding the optical nonlocality, while the perfor¬ 
mance of the deep-snbwavelength optical beam propagation with near-zero phase variation is 
also explored with respect to the variation of the material loss. In the following sections, the 
diffraction-free optical beam propagation with near-zero phase variation is discnssed in the¬ 
oretical analysis based on EMT and optical nonlocality, together with nnmerical simnlation 
in Sec. 2 , followed by the conclnsion in Sec. 3. 

II. THEORY AND DISCUSSION 

A. Effective medium theory and optical nonlocality 

Illnstrated in Fig. 1(a), the metal-dielectric mnltilayer stack is composite of alternating 
layers of silver (Ag) and titaninm pentoxide (TiaOs), with the thickness and a^, respec¬ 
tively. The permittivity of Ag follows the Drude model Sm = ^oo — ^p! + with 

the permittivity constant Soo = 5.7, the plasma freqnency Up = 1.37 x 10^® rad/s, and the 
damping factor 7 = A ■ (8.5 x 10^^) rad/s. According to the previons experimental demon¬ 
stration [35], the material loss of Ag can be rednced by the annealing process so that the 
damping factor 7 is decreased. Here the damping factor ratio A is included in order to tune 
the material loss of Ag and is set to be A = 0.1 in the following analysis. The permittivity 
of TiaOs is Ed = 5.83. When the thickness of each layer is sufficiently small, i.e., \ka\ -C 1, 
where k is the wave vector in the layer and a is the thickness of the layer [36], the multilayer 
stack can be regarded as a homogeneous anisotropic effective medium with a permittivity 
tensor based on the EMT 
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which is related to the filling ratio of Ag and TisOs, fm = am/ (am+Od) and fd = ad/ ( 0 ^+ 0 ^). 
According to Eq. (1), it is interesting that if the real part of the permittivity of Ag is 
greater than that of the imaginary part, i.e., Re(e:m) ^ Ini(£m )5 fhe EMT-based ENZ 
wavelength (where Re(£™"'") = 0) and the EMT-base ENP wavelength (where the maximum 
of Im(e™"'") is obtained) can be approximately associated with the EMT-based filling ratio 
of Ag as Ed/ied-Re{em)) and -Re(em)/(^d - Re(£m)), respectively. 

Therefore, the EMT-based ENZ-ENP wavelength is related to the approximate condition 
of Ed = —Re(£m); which implies that once the EMT-based ENZ wavelength overlaps with 
the EMT-based ENP wavelength, the EMT-based filling ratio of Ag should be Z™"*" ~ 1/2. 
Further algebraic calculation indicates that the accuracy result of the EMT-based filling 
ratio of Ag is = 0.500, with respect to the permittivity of Ag in terms of Drude 

model and the permittivity of TiaOs. Correspondingly, the EMT-based effective permittivity 
tensor reads e™"*" = 0.000 -|- 0.012i and = 7.238 -|- 2798.497i, presenting an extremely 
anisotropic permittivity property of the multilayer stack. Regarding a TM-polarized light 
with non-vanishing Hy, and field components propagating along the z-direction in 

the x-z plane, the iso-frequency contour (IFC) of the multilayer stack reads 


,2 / UM i , 7 2 / nivrr _ ,2 

'^x/^z '^z/^xy — 


( 2 ) 


while considering the extremely anisotropic permittivity property, i.e., S> 1 3> 

~ 0 at the EMT-based ENZ-ENP wavelength, the IFC equation in Eq. (2) can 
be approximately written as 

k-Jko = J e™’' (1 - (k^/kaf 

A- I - , P) 

Therefore, as shown in Fig. 1(b), according to Eq. (3), the IFC of the multilayer stack at 
the EMT-based ENZ-ENP wavelength forms fiat lines of the value kz/k^ ^ with 

respect to a wide k^/ko wave vector range, which leads to a diffraction-free optical beam 
propagation in the multilayer stack since all spatial components propagate with the same 
phase velocity along the ^-direction. Furthermore, due to the fact that Re(£™"'") = 0 and 
^ phase variation of the diffraction-free optical beam propagation in the 

multilayer stack is close to zero. 

However, previous studies show that the metal-dielectric multilayer stack possesses strong 
optical nonlocality [36-39], leading to the effective permittivity tensor not only related to the 
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frequency but also to the wave vector, especially when the frequency approaches to the ENZ 
position of the structure. Such optical nonlocality can be analyzed via the transfer-matrix 
method, in which the multilayer stack is considered as a one-dimensional photonic crystal 
with the dispersion relation to the TM-polarized light as 


cos(/C;j(am -h ad)) = cos{kmam) cos{kdad) 

i^^mkd/£dkm ^dkm/^mkd)/‘^ sin(/CffjOj^) sin(/C(^Oi^), 


( 4 ) 


with the wave vector ki = \/Sik^ — k^ for i = m,d. Accordingly, the effective permittivity 
component in x- and y-direction including the optical nonlocality reads 


^nonioc ^ arccos^[cos(v/^A;oam) cos{y/e^koad) 




sin(y^A;oam) sin(y^/coarf)]/(A;o(a„ -F a^)^). 


( 5 ) 


On the other hand, the nonlocal effective permittivity component in the 2 :-direction can be 
analytically obtained according to the effective permittivity dehnition = {Dz)/{Ez) 

based on the transfer-matrix method. Due to the optical nonlocality, the nonlocal ENZ 
wavelength and the nonlocal ENP wavelength for the multilayer stack will be different from 
the EMT-based ENZ-ENP wavelength. Therefore, modihcation on the multilayer geometry 
is necessary for a certain period of -|- ad, in order to exactly overlap the nonlocal ENZ 
wavelength with the nonlocal ENP wavelength. Figure 1(c) shows an example about the 
nonlocal effective permittivity components for a multilayer stack with the period of a^+ad = 
20 nm. Compared with the EMT results, the nonlocal results clearly indicate the shift of 
the ENZ-ENP wavelength from the EMT calculated position of A™"*" = 466.896 nm to the 
new position of A°°“*°'^ = 468.424 nm, while the nonlocal effective permittivity components 
are = 0.000 -|- 0.012i and = 7.214 -|- 2800.879i Furthermore, the hlling ratio of 

Ag related to the optical nonlocality (nonlocal hlling ratio for short) also slightly changes to 
= 0.501, with respect to the period of the multilayer stack. In more detail. Figs. 2(a- 
b) individually display the variations of the hlling ratio of Ag and the ENZ-ENP wavelength 
with respect to diherent periods of the multilayer stack caused by the optical nonlocality. 
For comparison, the EMT results are also plotted. It is clear that as the multilayer stack 
period gets larger, the nonlocal hlling ratio of Ag increases and the nonlocal ENZ-ENP 
wavelength also shifts to a longer wavelength. 
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B. Diffract ion-free and near-zero phase variation propagation 


Based on the finite element method (FEM), the diffraction-free optical beam propagation 
with near-zero phase variation is demonstrated according to the nnmerical simnlation in 
terms of two anti-phase TM-polarized Ganssian beams that are excited via the stimulating 
port boundary condition with an initial electric field distribution following the Gaussian 
distribution function as exp (—(x -|- cqY/wq), where cq determines the beam center and wq 
determines the waist size. Otherwise, the scattering boundary condition is applied to prevent 
the unwanted scattering from other boundaries of the simulation region. The simulation 
results are presented as the amplitude of the electric field, normalized by the maximum 
value. Figure 3 shows the calculated optical beam propagation of the Gaussian beams with 
ultra-narrow beam waists of wq = 25 nm in the homogeneous anisotropic effective medium 
at the EMT-based ENZ-ENP wavelength of A™"*" = 466.896 nm. Besides, Fig. 4 shows 
the optical beam propagation results calculated in the 20 nm-period multilayer stack at the 
nonlocal ENZ-ENP wavelength = 468.424nm. As shown in Figs. 3(a) and 4(a), it 

is clear that the deep-subwavelength Gaussian beams can propagate over a long distance 
(> 2A°°“*°'^) without any wave front distortion. At the same time, near-zero phase variation 
is also obtained along the propagation direction. The centre-to-centre distance of 120 nm for 
the two deep-subwavelength Gaussian beams remains well-defined as the beams propagating 
across the multilayer stack from the bottom to the top. It is noted that the confinement 
and propagation of the deep-subwavelength Gaussian beam inside the multilayer stack is 
entirely due to the unique extremely anisotropic permittivity property, and the optical beam 
path is well determined by the beam launching position, which is distinguished from the 
situation in a subwavelength optical waveguide [17], where the optical mode is confined 
by the waveguide boundary. Besides the straight optical beam propagation, the flow of 
light can be flexibly modeled through controlling the local geometry of the multilayer stack. 
For instance, the beam path can be manipulated by gradually varying the direction of 
multilayers, since the direction of the optical beam propagation is always vertical to the 
multilayer interface. The designed geometries for achieving 90° and 180° bending of deep- 
subwavelength Gaussian beams and the simulation results are shown in Figs. 3(b-c) and 
4(b-c), for the effective medium and the multilayer stack, respectively. It is noted that for 
the effective medium the anisotropic effective permittivity tensor is transformed according 
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to the geometric relation between the global x|/z-coordinates and the local nwrc-coordinates 
based on the transfermation optics as 

o, , (Eu cos'^ 9 + e^sin^ 9 sin 6^ cos 6^ \ ^ ^ 

^ = L , 2 9 • 

\ [Eu — Ew) sin 9 cos 9 Eu sin 9 + e^ cos 9 I 

The results indicate that the flow of light can indeed be manipulated while maintaining 
the diffraction-free propagation with near-zero phase variation. For comparison, Fig. 4(d) 
presents the optical beam propagation at the wavelength of A = 464 nm, which is off the 
nonlocal ENZ-ENP wavelength of the 20 nm-period multilayer stack, while the correspond¬ 
ing nonlocal effective permittivity components read = 0.121 -|- 0.012i and 

gnonioc| _ —335 266 -|- 42.895h It is clear that by overlapping the ENZ wavelength 

and the ENP wavelength, the multilayer stack with such an extremely anisotropic effective 
permittivity tensor supports diffraction-free propagation with near-zero phase variation. 
Also, it implies that the diffraction-free propagation supported by the anisotropic effective 
permittivity tensor is not extremely sensitive to the value of wavelength, but for the near-zero 
phase variation, the ENZ wavelength must be satisfied. 

In the previous analysis, the damping factor ratio of Ag is set as A = 0.1 in order to 
maintain a small material loss to extend the propagation distance of the Gaussian beams. 
In fact, the variation of the damping factor in the Drude model will affect the optical 
beam propagation loss and the phase variation. Eigure 5(a) represents the variation of the 
propagation length and the effective refractive index along the propagation direction with 
respect to different damping factor ratios of Ag. The effective refractive index is defined 
based on the nonlocal effective permittivity component as according to 

Eq. (5). The real part of the effective refractive index Re(n 2 ) is proportional to the phase 
variation of the propagating beam, while the imaginary part \m.{nz) is inversely proportional 
to the propagation length of the beam as A/(47rIm(n^)), where A is the wavelength of the 
beam in free space. It is clear that as the damping factor ratio of Ag gets larger, the real 
part of the effective refractive index increases, leading to an obvious beam phase variation 
along the propagation direction. On the other hand, the propagation length also decreases 
as the material loss gets higher. Figures 5(b-d) show the propagation of the two anti-phase 
Gaussian beams in the 20 nm-period multilayer stack with different damping factor ratio of 
A = 0.1, 0.5, and 1.0 at the nonlocal ENZ-ENP wavelengths, illustrating the dependences 
of propagation length and phase variation on the damping factor ratio of Ag. It is shown 
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that a smaller material loss will benefit the diffraction-free and near-zero phase variation 
optical beam propagation in the multilayer stack. 

Finally, it is worth mentioning that according to the previous experimental work of op¬ 
tical hyperlens [11], our above theoretical analysis of extremely anisotropic metamaterials 
can be realized in experiments. The ultra-narrow Gaussian beams can be achieved by in¬ 
scribing nanometer scale slots or holes into a thin chrome layer coated on the top of the 
sample surface, while the impedance matching substrate can be applied at the bottom of 
the multilayer stack in order to reduce the reflection caused by the mismatched impedance. 


III. CONCLUSIONS 

To conclude, extremely anisotropic metal-dielectric multilayer metamaterials with identi¬ 
cal ENZ wavelength and ENP wavelength have been designed based on the transfer-matrix 
method including the optical nonlocality. Such metamaterial possesses a near-flat IFC with 
wave vector values close to zero over a broad wave vector range. This unique property is 
utilized to obtain the diffraction-free deep-subwavelength optical beam propagation with 
near-zero phase variation, and it is demonstrated that the optical beam propagation can 
be manipulated flexibly by tuning the direction of the multilayer stack. The influences of 
multilayer period and material loss on the optical beam propagation are also studied. The 
current study is of great potentials in many applications such as optical imaging, optical 
integration, and on-chip optical communication. 
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FIGURE CAPTIONS 


FIG. 1. (Colour online) (a) Schematic of Ag-TiaOs multilayer stack, (b) The IFC of 
the Ag-TisOs multilayer stack at the EMT-based ENZ-ENP wavelength, where the wave 
vector kz/ko possesses near-zero flat real part (red lines) and imaginary part (black lines), 
compared with the IFC of air (green circle), (c) The EMT-based (dashed curves) and 
the nonlocal (solid curves) effective permittivity components, with the marked EMT-based 
(dashed vertical line) and the nonlocal (solid vertical line) ENZ-ENP wavelength. 

FIG. 2. (Colour online) The variation of (a) the Filing ratio of Ag and (b) the ENZ-ENP 
wavelength with respect to different periods of the Ag-TiaOs multilayer stack with optical 
nonlocality, compared with the EMT-based results representing the zero-period case. 

FIG. 3. (Colour online) The diffraction-free optical beam propagation with near-zero phase 
variation of two anti-phase deep-subwavelength TM-polarized Gaussian beams in the effec¬ 
tive medium for the case of (a) straight propagation, (b) 90° beam bending, and (c) 180° 
beam bending at the EMT-based ENZ-ENP wavelength. 

FIG. 4. (Colour online) The diffraction-free optical beam propagation with near-zero phase 
variation of two Gaussian beams in the Ag-TiaOs multilayer stack for the case of (a) straight 
propagation, (b) 90° beam bending, and (c) 180° beam bending at the nonlocal-based ENZ- 
ENP wavelength, (d) The optical beam propagation of two Gaussian beams off the nonlocal- 
based ENZ-ENP wavelength. 

FIG. 5. (Colour online) (a) The variation of the propagation length (blue curve) and 
the real part of the effective refractive index (red curve) with respect to different damping 
factor ratios of Ag. (b-d) The diffraction-free optical beam propagation with near-zero 
phase variation of two Gaussian beams in the Ag-TisOs multilayer stack with respect to the 
damping factor ratio of (b) A = 0.1, (c) A = 0.5, and (d) A = 1.0 at the nonlocal-based 
ENZ-ENP wavelengths. 
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